Signatures of polaronic charge ordering in optical and dc conductivity using 

dynamical mean field theory. 
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We apply dynamical mean field theory to study a prototypical model that describes charge order- 
ing in the presence of both electron-lattice interactions and intersite electrostatic repulsion between 
electrons. We calculate the optical and d.c. conductivity, and derive approximate formulas valid 
in the limiting electron-lattice coupling regimes. In the weak coupling regime, we recover the usual 
behavior of charge density waves, characterized by a transfer of spectral weight due to the opening 
of a gap in the excitation spectrum. In the opposite limit of very strong electron-lattice coupling, 
instead, the charge ordering transition is signaled by a global enhancement of the optical absorption, 
with no appreciable spectral weight transfer. Such behavior is related to the progressive suppression 
of thermally activated charge defects taking place below the critical temperature. At intermediate 
values of the coupling within the polaronic regime, a complex behavior is obtained where both 
mechanisms of transfer and enhancement of spectral weight coexist. 



I. INTRODUCTION 

As opposed to conventional charge density wavesji 
that are well understood in terms of lattice-driven in- 
stabilities of the Fermi surface in metallic systems, 
there is no unified description of the charge ordering 
(CO) transitions observed in strongly interacting sys- 
tems, or "bad metals". Typical examples of such sys- 
tems are transition-metal oxides, which are invariably 
characterized by a complex interplay between several 
microscopic interactions, involving the charge, spin, or- 
bital and lattice degrees of freedom. This complexity 
seems to preclude the identification of a simple, com- 
mon, charge ordering mechanism. Still, electron-lattice 
interactions are always present to some extent, and of- 
ten play a dominant role in driving the CO transition 
in these systems. Charge ordering phenomena involv- 
ing a strong electron-lattice coupling have been found 
in manganitesj2>2»a nickelates;^ layered cobaltates^ 
magnetite^ ' 10 ' 11 vanadates ; 12 ' 13 oxoborates, 1 - as well as 
in other inorganic low-dimensional systems i 15 ' 16 The 
lattice degrees of freedom could also be relevant in 
the charge ordered phases of the so-called "telephone- 
number" ladder compounds^ as well as in two- 
dimensional organic salts, 18 although in those cases most 
theoretical interpretations up to now have focused on 
models with purely electronic interactions. 

In this work we focus on a minimal model which de- 
scribes charge ordering in the presence of electron-lattice 
interactions, with particular attention to the polaronic 
regime obtained at strong coupling. The model consists 
of electrons on a bipartite lattice, at a commensurate con- 
centration of one electron on every two sites, interacting 
locally with dispersionless lattice vibrations. The elec- 
trons also interact mutually via an intersite electrostatic 



repulsion, which can be thought of as the screened part 
of the long-ranged Coulomb potential. While this over- 
simplified model only retains part of the complex physics 
involved in real systems, its solution can be helpful to 
clarify certain aspects of the charge ordering phenomena 
that are common to systems with strong electron-lattice 
interactions, and that can in principle be identified ex- 
perimentally. To be specific, we solve the model by per- 
forming the following approximations: 

(i) The magnetic degrees of freedom are taken out of 
the game by resorting to spinless electrons. This enforces 
locally the constraint of no double occupancy character- 
istic of the limit of strong "Hubbard" (on-site) repulsion. 
It is appropriate for our purposes, as long as we do not 
aim at describing the effects of electronic correlations on 
the low-energy physics, such as the existence of a quasi- 
particle peak or the induced magnetic exchange. In prin- 
ciple, this approximation is viable as long as the magnetic 
and charge degrees of freedom are governed by different 
energy scales. In such case, the magnetic interactions 
are not expected to have a strong influence on the charge 
ordering pattern (although exceptions do exis t 19 ' 20 ' 21 ). 
This separation of energy scales is realized for instance 
in magnetite (FesO^r^- where the CO transition occurs 
within a ferromagnetic phase that preexists at a much 
higher temperature. More examples can be found in sys- 
tems (such as VsOgi 3 -) where the CO transition takes 
place within a paramagnetic phase, the magnetic order 
setting in at a much lower temperature. The spinless 
approximation is also relevant to half-metals such as the 
colossal magnetoresistance manganitcs. 

(ii) We consider a single electronic band, which rules 
out the orbital degrees of freedom and their possible or- 
dering, that is known to play an important role in specific 
compounds. 
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(iii) The lattice vibrations are assumed adiabatic, 
i.e. their characteristic energy is much smaller than 
the electronic bandwidth, which is typically the case in 
transition-metal oxides. The adiabatic approximation is 
enforced here by treating the lattice degrees of freedom 
as static variables with a given (thermal) statistical dis- 
tribution. As a consequence, the low-energy spectral fea- 
tures that derive from the quantum nature of the phonons 
are lost. These would appear at energies comparable to, 
or below the phonon energies. Such features are any- 
how of minor quantitative importance in the polaronic 
regime of interest here, where most of the spectral weight 
is moved to higher energies, of the order of the polaron 
binding energy. The same argument justifies the neglect 
of low-energy phenomena originating from electronic cor- 
relations as discussed at point (i) above, in all cases where 
the polarons set the dominant energy scale. 

(iv) We apply single-site Dynamical Mean Field The- 
ory (DMFT), appropriately adaptedi&^^i^ to account 
for the charge unbalance between neighboring sites in the 
charge ordered state. This approach is known to deal 
very effectively with local interaction mechanisms, re- 
gardless of their strength. It is therefore expected to give 
an accurate description of the electron-lattice physics. As 
was shown in Ref. l23l this approach correctly describes 
the crossover between the weak coupling charge density 
wave regime and strong-coupling polaronic regime, and 
allows to shed light on the role played by defects of the 
lattice polarization at the ordering transition. A lesser 
accuracy is achieved in treating the inter-site Coulomb re- 
pulsion, since non-local interactions reduce to the mean- 
field (Hartree) level in single-site DMFT. This is how- 
ever not a primary limitatio n 19 ' 20 in the regime we are 
mainly interested in, where the physics is dominated by 
the electron-lattice interaction. 

The main result of this work is the identification of 
qualitatively different behaviors of the electrodynamic 
response in the two fundamentally distinct regimes of 
charge ordering, governed by the strength of the electron- 
lattice interaction. In the weak coupling regime, the elec- 
trodynamic response is dominated by the opening of an 
energy gap in the excitation spectrum of an otherwise 
metallic system. At the ordering transition, the optical 
spectral weight is transferred from the Drude peak to 
the region above the gap, and the d.c. conductivity ac- 
quires a semiconducting character. On the contrary, in 
the strong coupling limit, charge localization is already 
present in the normal phase, due to the formation of self- 
localized polarons. Correspondingly, the low energy spec- 
trum is already strongly depleted by a polaronic "pseudo- 
gap" , and the genuine gap opening occurring at the tran- 
sition does not give rise to any appreciable depletion or 
shift of the low frequency spectral weight. Instead, the 
ordering of polarons results in a global enhancement of 
the optical spectral weight, which is opposite to what is 
expected from a conventional charge density wave. The 
two scenarios described above coexist in the polaronic 
regime at moderate values of the electron-lattice cou- 



pling. This leads to a complex behavior presenting both 
a transfer and an enhancement of the spectral weight, 
resulting in a non-monotonic temperature dependence of 
the total optical weight within the CO phase. The in- 
clusion of a direct intersite Coulomb interaction extends 
this complex intermediate region to larger electron-lattice 
interaction strengths. 

The solution for the one-particle spectral function has 
already been published in Ref. |23| . in a model with 
electron-lattice interactions alone. Here we extend that 
treatment to calculate the optical and d.c. conductiv- 
ity throughout the parameter space, including an inter- 
site Coulomb repulsion between electrons. We also de- 
rive analytical expressions for these quantities that are 
valid in the polaronic regime, elucidating the role played 
by charge defects in both optical absorption and charge 
transport. It can be noted that an analogous model, 
without the Coulomb interaction term, was studied vari- 
ationally in Ref. H|| There, however, the physics of de- 
fects was not addressed. 

The paper is organized as follows. The general for- 
malism for treating the charge ordered, broken-symmetry 
phase within the DMFT is presented in Sec. II. The equa- 
tions are solved numerically in Sec. III. The optical and 
d.c. conductivity are calculated in Sec. IV and V respec- 
tively, making use of the Kubo formula. The results are 
discussed in Sec. VI. 



II. MODEL AND SELF-CONSISTENT 
SOLUTION 

We consider the following Holstcin-Coulomb Hamilto- 
nian: 

H = -±J2(ctc J + h.c.) + ±kX? (1) 

v Z <ij> 

i <ij> 

where the operator cf(ci) creates (destroys) a spinless 
electron at site i, n is the average electron density per site 
and Xi are the displacements of local oscillator modes, 
that are coupled to the local electron density via the pa- 
rameter g. In addition to this electron-phonon interac- 
tion a la Holstein, the electrons interact mutually via a 
nearest neighbor repulsion term V of Coulomb origin. 
The scaling of the repulsive term is chosen in such a way 
that the energy cost for a deviation from perfect charge 
ordering (i.e. flipping an occupied and an empty site) 
equals V. The scaling of the kinetic term t/y/z yields 
a finite free-electron bandwidth in the limit of infinite 
connectivity (number of nearest neighbors z — > oo). We 
shall consider for simplicity a semi-elliptical density of 
states (DOS) of half-width D — 2t, corresponding to a 
Bethe lattice of infinite connectivity. Nevertheless, the 
present calculation scheme can be easily generalized to 
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other model DOS, to include more realistic band struc- 
tures as obtained for example from ab-initio calculations 
of specific systems. We shall specialize to the commen- 
surate concentration n = 1/2, which is equivalent to a 
quarter-filled band in the case of electrons with spin. In 
the present spinless case, the chemical potential is fixed 
at u = from particle-hole symmetry. 

An order parameter for the CO transition can be de- 
fined as the charge disproportionation 



An = 7ia — Ub > 



(2) 



which varies between and 1, the average density being 
given by n — (ua + ub)/2 = 1/2 (here A and B label re- 
spectively the charge- rich and charge-poor sublattices). 
The formalism needed to treat the electron-phonon in- 
teraction was set up in Ref. [23| by mapping the origi- 
nal lattice problem onto a pair of coupled impurity mod- 
els. In the limit of infinite connectivity, any additional 
non-local interaction term self-averages to its mean-field 
(Hartree) value. Therefore, the extra Coulomb term can 
be straightforwardly included via a site dependent chem- 
ical potential shift in the Weiss fields, i. e. the local 
fields that take into account the effects of the electron 
itinerancy. On the Bethe lattice, the new Weiss fields 
read: 



(G- 1 ) 



(Go 1 )" = 



- —An - —G BB 
2 4 

—An - —G AA . 
2 4 



(3) 
(4) 



It can be observed by direct inspection of Eqs. (J3]|l|) that 
the Coulomb interaction does not affect the properties of 
the system in the normal phase (where An = 0), which 
is a drawback of the present approximation. 

The above system of equations is closed self- 
consistcntly by expressing the site-diagonal propagators 
G AA and G BB in terms of the self-energy arising from 
the local interaction with the static phonon field: 



G 



AA 



dX 



Pa(X) 



(G^) A -gX : 



(•5) 



and similarly for G BB . Here P A (X) and P B (X) are the 
probability distribution functions (PDF) for the phonon 
displacements on the A and B sublattices, which can be 
obtained by tracing out the electronic degrees of freedom 
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p A (x)oce^ kx2 ^Y[[(G- 1 ) A - 9 x} e* 



(6) 



The numerical solution is achieved by successive itera- 
tions of the above equations, starting from a given ansatz 
for the sublattice propagators G AA and G BB . In the fol- 
lowing we shall take the half-bandwidth D as the unit of 
energy. 

We now briefly analyze the phase diagram obtained 
from the DMFT solution of the model Eq. {T). Fig. [1] 




FIG. 1: Critical temperature delimiting the CO phase as 
a function of A = g 2 /2kD, for different values of V/D. 
Shaded areas below T c mark polaronic regions where de- 
fects are present. The existence of defects is tracked by 
looking at metastable minima of the adiabatic potentials 
Va,b = -riog(P A , B ). 



shows the critical temperature for the charge ordering 
transition as a function of the electron-phonon coupling 
parameter A = g 2 /2kD, defined as the polaron energy 
Ep = g 2 /2k in units of the half-bandwidth. At V — 0, we 
recover the bell-shaped result obtained in Ref. [23]. Such 
bell shape originates from the different mechanisms that 
drive the charge ordering in the two limiting regimes of 
small and large A. At weak coupling, the critical temper- 
ature monotonically increases with the electron-phonon 
coupling, i.e. the same qualitative trend predicted by the 
usual BCS theory^ In this regime, the charge dispropor- 
tionation and the long-range order take place simultane- 
ously at T c . At strong coupling instead, due to polaron 
formation, a local charge segregation corresponding to the 
localization of the electrons on randomly distributed sites 
occurs at a temperature T w Ep , much larger than order- 
ing temperature T c itself. The actual critical temperature 
T c marks the onset of spatial ordering of such randomly 
distributed polarons. It is proportional to the "charge 
superexchange" J = D / (4A) and therefore decreases as 
the electron-phonon coupling increases. The asymptotic 
strong coupling expression for the critical temperature is 
T c = D/ (16A). From the above discussion we see that the 
maximum of T c , obtained at A ~ 1 for V = 0, signals the 
crossover between the weak-coupling and the polaronic 
behavior. 

The ordered phase below T c is also qualitatively differ- 
ent depending on the value of the electron-phonon cou- 
pling. In particular, in the strong coupling regime, the 
preexisting polarons progressively order upon lowering 
the temperature, and only attain a perfect alternate or- 
der at T = 0. At any finite temperature, a certain num- 
ber of defects exist within the CO phase, in the form 
of charges localized on the "wrong" sublattice. The ex- 
istence of such thermally-induced defects is intimately 



4 



2 
1.8 
1.6 
1.4 
1.2 

1 

0.8 
0.6 
0.4 
0.2 




-3 -2 



V=0 t/T>0.3 
V=0 T/T.=0.9 
V=1 T/T>0.9 



FIG. 2: Sublattice phonon PDF at A = 2. The different curves 
are for V = (T/T c = 0.3 and 0.9) and V = 1 (T/T c = 0.9). 
Points correspond to the DMFT self-consistent solution, lines 
represent the strong coupling theory Eq. (|B5() . The locus and 
shape of the peaks does not depend explicitly on V . 




FIG. 3: Evolution of the charge disproportionation for A 
0.4 (a) and for A = 2.0 (b), at different values of V. 



related to the existence of local lattice distortions in the 
normal phase, i.e. polaron formation^ We can therefore 
identify the existence of a polaronic CO phase, as opposed 
to the weak-coupling charge density wave, as the phase 
in which these defects are present. 

The defect phase can be quantitatively characterized 
by analyzing the sublattice PDF's, which acquire a bi- 
modal structure for sufficiently large A. More efficiently, 
one can look for the existence of two non-degenerate 
minima in the sublattice adiabatic potentials denned as 
Va,b = — T log Pa, b- The result of this procedure is il- 
lustrated by the shaded areas in the phase diagram of 
Fig. [U We see that the defect phase emerges right at the 
polaron crossover, i.e. in the region below the maximum 
of T c vs A. 

Including a direct intersite repulsion clearly favors both 
types of charge ordering, as it increases T c for all values 
of A (compare the different curves in Fig. [1]). Still, the 
evolution of the phase diagram with V suggests a non- 
trivial interplay between the direct Coulomb repulsion 
and the polaronic physics, as a finite V shifts the polaron 
crossover (as well as the boundary of the defect phase) 
towards larger values of A. It can be argued that, by 
treating the intersite Coulomb interaction beyond mean- 
field, charge fluctuations of the same nature of the lattice 
defects evidenced above could start playing a role^ For 
this reason, in the following we shall restrict to the regime 
of small to moderate V, where the dominant physics is set 
by the electron-phonon interactions, therefore suppress- 
ing the fluctuations of the Coulomb interaction term. 

The sublattice PDF P A (X) evaluated at A = 2, well 
inside the polaronic regime, is shown in Fig. [21 Defects 
in the A sublattice give rise to a minority peak of op- 
posite "polarization" with respect to the preferred one. 
The area of the minority peak can be used to define the 
number of defects (rid) in the CO phase. In the strong 
coupling regime, it tends to 1/2 at T = T c , where the 
two peaks become equivalent (polarons are randomly dis- 



tributed over the two sublattices), while it vanishes ex- 
ponentially at T = 0, when perfect charge ordering is 
achieved. This occurs because, as the temperature is re- 
duced, the population of the metastable minimum of the 
potentials Va,b is progressively depleted, leading to an 
exponential reduction of n<j. It should be noted that the 
possibility of strictly rid = is related to our classical 
treatment of the lattice degrees of freedom, which breaks 
down at temperatures much smaller than the characteris- 
tic phonon energies. Properly including phonon quantum 
fluctuations 2 ^, would lead to the appearance of charge de- 
fects even in the zero temperature limit. 

In the limit A — > oo, the number of defects is related 
to the order parameter through 

An = 1 - 2n d . (7) 

At finite A in the polaronic phase (A > 1), the above Eq. 
© has to be replaced by An = (l-2n d )[l-l/(8A 2 ) + . . .], 
which properly accounts for the fact that the charge dis- 
proportionation An does not strictly tend to 1 even when 
the polarons are perfectly ordered (rid = 0). This is be- 
cause at finite A the electronic wavefunction is not fully 
localized and acquires a finite extension on the neigh- 
boring sites, which necessarily pertain to the minority 
sublattice. 

The evolution of the charge disproportionation with 
T in different regions of the phase diagram is reported 
in Fig. [3] In the polaronic regime, An is a universal 
function of T/T c , independent of V, as shown in Fig. 
[3]b. This means in particular that the number of de- 
fects rid is a function of T/T c alone, and V only en- 
ters implicitly through the determination of T c . Inter- 
estingly, the gap ratio At=o/T c can become very large 
in the polaronic regime, due to the preexisting pseudo- 
gap A cx Ep ^ T c . 23 In the opposite weak electron- 
phonon coupling limit, we recover a BCS-like behav- 
ior for An. In that case the gap is given by A = 
(XD + V/2)An. Its value at T = satisfies the mean-field 
ratio At=q/T c = 3.54. It is worth noting that this limit 
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is only reached asymptotically as A — > 0. For small but 
finite A, T c is dramatically reduced by thermally induced 
lattice fluctuations , 23 ' 30 while At=q is not. Therefore the 
gap/T c ratio rapidly increases with A even well outside 
the polaronic region. 



III. OPTICAL CONDUCTIVITY 

A. DMFT formulation 

We now derive the Kubo formula for the optical con- 
ductivity in the CO phase. To this aim we develop a gen- 
eral formalism which makes explicit use of translational 
invariance as in the hypercubic lattice, following the lines 
of Ref. HJ. In this framework it is possible to write ex- 
pressions for the electron propagators in fc-space. This is 
achieved through a canonical transformation that defines 
new electron creation operators c^' B — (ck ± Ck+Q)/\/2, 
where Q = (jt, it, •••) is the instability wave vector in 
any dimensions and k spans the reduced Brillouin zone 
(RBZ). With the above transformation, the tight bind- 
ing term becomes J2 k £k(c k A c k J + c\fc£), where £k is 
the original non-interacting band dispersion. The fully 
interacting fermion propagators on the bipartite lattice 
can be defined as the matrix elements 



G 



af3 



-«<Tc£(i)c£(0)) (a,0 = A,B) 



(8) 



of the 2x2 matrix G k . The corresponding spectral func- 
tions are given by p(k, to) = —ImG^uj)/?:. Once the 
local self-energies on the two sublattices are known from 
the solution outlined in Sec. [ill t ne Green's functions of 
Eq. ([5]) are obtained by inverting the matrix 



Gk 1 



(9) 



with z A = uj + i8 + p — Y,a(uj), and similarly for z B . 

In the tight binding model of Eq. ([1}, the current 
operator along a given (say x) direction reads 



jr.. I , 



EM+^ 



(10) 



where the sum extends over all sites i of the lattice and 
their z nearest neighbors, identified by the vectors 5. 
Transforming to the sublattice operators, the current op- 
erator can be expressed as: 



In a hypercubic lattice, only the two neighbors along the 
x direction contribute to the above sum, and the corre- 
sponding current vertex is «k 



2^j sinfc x . 



In the context of DMFT, due to the vanishing of vertex 
correction a 32 ' 33 ' 34 the current-current correlation func- 
tion can be expressed exactly in terms of single particle 
Green's functions. This simplification still holds in the 
broken-symmetry phase because the k — > — k symme- 
try is preserved even in the reduced Brillouin zone. Re- 
minding the definitions Eq.®, we can write the current- 
current correlation function as 



(TJ x {t).J x {0))= E 



vl Tr 



T x G k (t)r x G k (-t) (12) 



keB,BZ 



where we have made use of the Pauli matrix t x to obtain 
a compact expression. The trace is performed in the 
sublattice indices A, B. 

The optical conductivity follows from the Kubo for- 
mula, upon transforming the current-current correlation 
function to the frequency domain: 



a{uS) =<7o deN(e)<j)(e) x 

J RBZ 



dv Tr [T x p(e, lo + v)T x p(e, v)] 



(13) 

f{y)-f{ V + U)) 



In the above equation, the constant o"o = Tre 2 a 2 /hv 
carries the dimensions of conductivity, a being the lat- 
tice spacing, v the volume of the unit cell, f(v) — 
[1 + e^" - ^)] -1 the Fermi function and p the spectral 
functions associated to sublattice propagators G. Tak- 
ing advantage of the local nature (fc-independence) of 
the self-energy, the sum over momenta in Eq. (fT2")) has 
been replaced by an integration over energies, weighted 
by the DOS N(e) and the current vertex </>(e) of the non- 
interacting lattice. Use has also been made of the fact 
that is invariant under the transformation k — > k + Q. 
The integration over the band dispersion e in the current- 
current correlation function can be performed analyti- 
cally, leading to: 



a(uj) = cr 



di,B(u + Ul v) K v + u )-fM. (14) 



J x 



keRBZ 



Ck 1 )- 



(ii) 



The function B is defined as: 



B(y + u>, v) 



4tt 



Re { X {z?, zf ; zt zi) - X (zt *?; zf, *)} 



(15) 
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where 



SI _ S2 I S2 SI 



and 



"2 

Si = 

6 = 



1/ • 



- z/ + i(5 — T, a (uj + v) 
id — Yl a [y) 



(17) 
(18) 

(19) 
(20) 



with a — A, B the sublattice index. The function JC is 
the Hilbert transform of the product N(e)<p(e), that can 
be evaluated analytically in the case of a semi-elliptical 
DOS. For this we take the following form of the current 
vertex ci(e) = (D 2 — e 2 )/3. This is chosen in such a way 
as to fulfill the f-sum rul o 35 ' 36 i 37 which relates the total 
optical spectral weight 



W 



duja(uj) 



(21) 



to the kinetic energy K of the interacting system. The 
demonstration of the f-sum rule in the broken-symmetry 
phase is explicitly carried out in Appendix C. 

In the following Sections we shall analyze separately 
the results in the weak electron-phonon coupling regime 
A < 1) m the extreme strong coupling limit A > 1, and 
in the most interesting polaronic regime A > 1. 



B. Weak coupling regime, A < 1 

In figure H] we report the results for the optical con- 
ductivity at A = 0.4, for different values of the intersite 
Coulomb interaction V. At such a moderate value of the 
electron-lattice coupling, although the overall behavior 
qualitatively agrees with what is expected in the con- 
ventional mean-field scenario (see Appendix A, and the 
dashed lines in Fig. [4|), the spectral features are apprecia- 
bly smoothened by the presence of lattice fluctuations. A 
measure of such broadening is provided by the variance 
s = \J2EpT of the energy fluctuations of the phonon 
field , 38 i 39 which is of the order of 0.1D in the example of 
Fig. HI This is a sizable fraction of the electronic band- 
width, which is however sufficiently small that the overall 
band picture remains qualitatively valid at temperatures 
comparable or below T c , which are the ones of interest 
here. 

Due to the classical approximation for the phonons, 
the Drude peak in the normal phase is replaced by an in- 
choerent gaussian peak around u> = 0^ whose width is 
proportional to s. Below the critical temperature, a finite 
charge disproportionation develops (cf. Fig. [3]) which 
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FIG. 4: (color online) Optical conductivity at A = 0.4, V — 
(a) and V = Q.2D (b). The curves in panel (a) are at equally 
spaced temperatures in the range T/T c = 0.4 — 2.9 (T c = 
0.026£>) and in panel (b) T/T c = 0.2 - 1.5 (T c = 0.052D). 
Curves corresponding to the disordered phase at T > T c are 
in red. The dashed line is the mean-field result Eq. ()A2|) at 
the lowest temperature. 



is reflected in the optical spectra through a progressive 
gap opening. Correspondingly, the low-frequency spec- 
tral weight is transferred to frequencies above the optical 
gap located at A opt ~ (2AZ? + V) An. The sharp square- 
root divergence at uj = A opt expected from the standard 
mean-field treatment (dashed line, corresponding to Eq. 
IA2|) is also smoothened due to the presence of thermal 
lattice fluctuations. These are also responsible for the 
subgap absorption tail observed in Fig. [4l which is ab- 
sent in the mean- field result. As shown in Ref. l26l . a 
similar broadening can also arise at T — due to quan- 
tum lattice fluctuations. 

The effect of a finite V is to strengthen the CO phase, 
without modifying the broadening of the spectral fea- 
tures. The optical gap increases both due to the in- 
creased value of the charge disproportionation An (see 
Fig©, and to the explicit contribution VAn, which rep- 
resents the extra electrostatic cost to move one particle 
from one sublattice to the other. 



C. Strong coupling limit, A 3> 1 

The results for the optical absorption in the strong cou- 
pling limit are shown in Fig. [S]for A = 2. The behavior 
is very different from the weak coupling regime analyzed 
previously. This is seen already in the normal phase, 
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FIG. 5: (color online) Optical conductivity at A = 2, V = 
(a) and V — 1 (b). The curves in panel (a) are at equally 
spaced temperatures in the range T/T c = 0.3 — 2.5 (T c = 
0.031D) and in panel (b) T/T c = 0.04 - 1.1 (T c = 0.2613). 
Curves corresponding to the disordered phase at T > T c are 
in red. The dashed line is the strong coupling result Eq. (|B7|) 
at the lowest temperature. 



where the Drude peak is replaced by a broad peak at fi- 
nite frequency. This peak arises due to the formation of 
small polarons, and reflects the optical transitions within 
the polaron internal structure. The peak position scales 
with the polaron energy roughly as w ~ 2Ep, and its 
width is again proportional to the spread s of the lattice 
fluctuations. 

Because the low-frequency spectral weight is already 
strongly suppressed in the normal phase, no clear gap 
opening is visible at T c . Indeed, in Fig. \5[a. (A = 2, 
V = 0) there is no visible depletion of the optical conduc- 
tivity at low frequency (apart from the natural evolution 
of the peak width governed by s) and no clear shift of 
the position of the polaron peak as the temperature is 
reduced below T c . Rather, the ordering transition leads 
to a sharp increase of the spectral weight associated to 
the polaronic peak. As we demonstrate hereafter, this in- 
crease is a direct consequence of the suppression of charge 
defects, and constitutes a distinctive signature of the po- 
laronic charge ordering. 

The observed behavior can be understood by noting 
that in the normal phase, only half of the sites neigh- 
boring a given polaron are unoccupied, and therefore 
available for an optical transition as induced by the cur- 
rent operator Eq. (| 10[) . The number of available empty 
neighbors increases as charge defects are progressively 
removed below T c , and so does the weight of the pola- 
ronic absorption peak, until each polaron becomes exclu- 
sively surrounded by unoccupied sites at T = 0. The 
above analysis can be carried out in general to show that 
the probability for polaronic optical transitions is propor- 
tional to UaO- — ub) (for polarons in sublattice A) and 
ns(l — Ua) (for polarons in sublattice B). The sum scales 



as (1 + An 2 ), leading to an increase of up to a factor of 
2 of the optical spectral weight in the CO phase. 

A more detailed discussion of the spectral weight en- 
hancement is provided in Appendix B, where we derive 
analytical expressions for the optical conductivity valid 
at large A. We report here the formula appropriate in the 
CO phase at V = 0, which is simply expressed in terms 
of the normal state polaron absorption a norm (uj) as 

er(cj) = 2[1 - 2n d (l - n d )]a norm (uj). (22) 

This equation directly relates the observed enhancement 
of the polaron peak with the suppression of the number 
of defects, in agreement with the qualitative argument 
given above, as can be easily shown using Eq. ([7]) above. 
Both this expression and its analog Eq. (|B8[) valid at 
V 7^ describe quite accurately the full DMFT results 
at A = 2, as illustrated in Fig. 

D. Spectral weight analysis 

The increase of spectral weight associated to the or- 
dering of polarons is best visualized by analyzing the 
integrated optical conductivity W defined in Eq. (|2"Tj) . 
The evolution of W vs. temperature is reported in Fig. 
O and exhibits a markedly different behavior at weak 
and strong electron-phonon coupling. In Fig. [SJc, cor- 
responding to A = 2 and V — 0, we see that upon low- 
ering the temperature W first decreases (in the normal 
phase) and then exhibits a sharp increase (in the ordered 
phase), saturating at T -C T c to a value which is about 
twice the normal phase value. Therefore, at large A the 
critical temperature corresponds to a minimum of W. 
This remains true in all the polaronic regime, down to 
the polaron crossover at A ~ 1. This is opposite to the 
usual behavior of charge density waves, where upon low- 
ering the temperature the optical weight first increases 
and then decreases, reaching a maximum at the critical 
temperature. Such conventional behavior is recovered at 
small A, as illustrated in Fig. [SJa for A = 0.4. 

These opposite behaviors can be understood in virtue 
of the optical sum rule demonstrated in Appendix C, 
which relates the total spectral weight to (minus) the to- 
tal kinetic energy K of the system. The latter is usually 
expressed as K = X^ 6 *-' 71 ^' wnere n e k is the momen- 
tum distribution function of the interacting electrons. As 
shown in Appendix C, an equivalent expression holds true 
even in the broken-symmetry phase [see Eq. (|C7[) ] , if one 
properly introduces the occupation number for the gen- 
eralized propagator G = 1/\^/zaZb — Plots of this 
quantity are shown in Fig [7] 

As pointed out in Ref. |4l|, any broadening of the step 
feature in n Ck around the Fermi energy, by populating 
states with higher energy, leads to an increase of K, and 
therefore to a decrease of the spectral weight W. In the 
case of a conventional charge density wave, illustrated in 
FigjTJa, the distribution n eh broadens both upon increas- 
ing the temperature above T c and upon opening a gap 
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FIG. 6: Total optical weight W as defined in Eq (|2TJ), for 
A = 0.4 (a), A = 1 (b) and A = 2 (c), at different values of V 
(the legend is indicated in panel b). Lines in panels (a) and 
(b) are guides to the eye, while in panel (c) they correspond 
to the analytical strong coupling approximation obtained by 
integrating Eqs. ()B7|I and ()B8|I down to a cutoff frequency 
equal to T. 



below T c . As a result, the critical temperature identifies 
a minimum of K , i.e. a maximum of W . 

The situation is different in the strong electron-phonon 
coupling regime (Fig. [7]b and c), where the formation 
of small (local) polarons involves all the states in the 
Brillouin zone. Correspondingly, the distribution n tk be- 
comes extremely flat. In this case, increasing the tem- 
perature above T c progressively destroys the local cor- 
relations that build up the polaron, and the momen- 
tum distribution eventually recovers some structure at 
the Fermi level, leading to a decrease of K (the effect 
is small because T c in this regime is much smaller than 
the energy scale Ep that governs the polaron dissocia- 
tion). The polaron ordering below T c also decreases K, 
because the emergence of two distinct sublattices effec- 
tively restores some momentum-space structure. In the 
polaronic regime, the critical temperature therefore de- 
fines a maximum of K, i.e. a minimum of W . 



FIG. 7: Momentum distribution function vs. band energy 
e = efc at temperatures at, below and above T c . The coupling 
strengths are the same as in Fig|6] and we have set V = 0. 

E. Polaronic regime, A > 1 

By following the evolution of the spectral weight from 
Fig. [S] a to c, it is interesting to see that at intermedi- 
ate values of the electron-lattice coupling (but still in the 
polaronic regime), the two competing trends evidenced 
above lead to a non-monotonic temperature dependence 
of the kinetic energy (and of the spectral weight) within 
the CO phase. This can be seen in particular the curve 
labeled V = at A = 1 in Fig. Ob. Note that a non- 
monotonic behavior of the kinetic energy analogous to 
the one described here was recently observed in the or- 
dered phase of the Falicov-Kimball model at intermediate 
UM. 

The spectral weight analysis presented above suggests 
that both features characteristic of the weak- and strong- 
coupling limits coexist in the polaronic regime at inter- 
mediate coupling strengths A > 1. We can check in Fig. 
[SJa (A = 1, V = 0) that such coexistence also natu- 
rally manifests in the optical spectra. Here, as in the 
strong-coupling limit of Fig. [5la, the polaronic nature 
of the carriers is testified by the presence of a broad 
finite-frequency peak already in the normal phase. In 
this case, however, where the polaron energy Ep is com- 
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FIG. 8: (color online) Optical conductivity at A = 1, V = 
(a) and V = 1.0 (b). The curves in panel (a) are at equally 
spaced temperatures in the range T/T c = 0.2 — 1.5 (T c = 
0.052D) and in panel (b) T/T c = 0.04 - 1.1 (T c = 0.25D). 
Curves corresponding to the disordered phase at T > T c are 
in red. The dashed line is the strong coupling result Eq. (|B7|) 
at the lowest temperature. 



parable with the free electronic bandwidth D, the peak 
position lies well below the strong-coupling estimate 2Ep. 
Also the shape and width of the peak deviate from the 
strong-coupling estimate, being both controlled by the 
ratio s/D (the usual symmetric gaussian shape is recov- 
ered as s ^> D, see Ref. |37l ). Most interesting however is 
the evolution of the optical absorption below T c , which 
shares similarities with both behaviors shown in Figs. UJa 
and a. Indeed, at intermediate values of the electron- 
phonon coupling, the optical absorption in the ordered 
phase exhibits both a marked enhancement and an ap- 
preciable transfer of spectral weight to higher frequencies. 
This composite behavior gives rise to an "isosbestic" re- 
gion at frequencies below the polaron peak, where the 
optical absorption is almost independent of temperature. 

It was observed by analyzing the phase diagram in Sec. 
II that the intersite repulsion V effectively pushes the sys- 
tem towards the weak electron-lattice coupling regime. 
This conclusion is also supported by the behavior of the 
optical conductivity spectra. For example, we see from 
Fig. [I] that at A = 1, a repulsion V = 1 is sufficient to 
move the system outside the polaronic phase. Accord- 
ingly, the polaron peak that was present at high temper- 
ature in Fig. [Sja has disappeared in the V = 1 spectrum 
of Fig. [Hb, and a more conventional gap opening is re- 
stored below T c . If instead one starts from the strong 
coupling value A = 2, an appreciable transfer spectral 
of spectral weight below the transition is recovered for 
V = 1, as well as a non monotonic behavior of the total 
spectral weight W (see Fig. [S]b and FigJSJc). 
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FIG. 9: d.c. conductivity for A = 0.4 (a), A = 1 (b) and 
A = 2 (c), for different values of V. Lines are guides to the 
eye, except for V — 1 in panels (b) and (c), where the line is 
the strong coupling formula Eq. ((24}. The full black line is 
the DMFT result obtained in the normal phase, and the thin 
dashed line is the standard polaronic expression Eq. (|23l> . 



IV. D.C. CONDUCTIVITY 

In this section we briefly describe the d.c. conductiv- 
ity, which is obtained from the previously derived Kubo 
formula Eq. (|13[) . by taking the limit u> — > 0. The numer- 
ical results for cr^.c. across the CO transition are shown 
in Fig. [9] for different values of V, at A = 0.4 (a), A = 1 
(b) and A = 2 (c) . These values are the same as analyzed 
in the preceding Section, and correspond respectively to 
the weak, intermediate and to the strong electron-lattice 
coupling regime. The CO transition can be clearly iden- 
tified as a knee in all the conductivity curves, except for 
A = 2, V — (in that case T c falls outside the range of 
Fig. [5Jc) In the weak coupling regime (Fig. [Ua), the 
opening of a gap changes the mobility from metallic-like 
to insulating-like, i.e. the slope of (Jd.c.{T) changes sign. 
The changes in the transport properties at the CO tran- 
sition become less marked in the polaronic regime. In 
this case the mobility is already thermally activated in 
the normal phase, and the CO transition only causes an 
increase of the activation barrier, which is reflected in a 
change of slope in the Arrhenius plots of Figs. [H]b and 
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Fig. lc. 

An expression for the d.c. conductivity valid in the 
strong coupling limit is derived in Appendix B. In the 
normal phase, one recovers the usual formula for the po- 
laronic activated behavior: 

„norm(rp\ "0 „ — Ep/2T /OQA 

° d - {T) = 32VTOT3/ 2 e ' (23) 

At temperatures below T c , this generalizes to 

* d . c .(T)=4n d (l-n d )aZ rm (n (24) 

We see that the ordering transition reduces the con- 
ductivity in a way that is once again related to the 
suppression of defects, being proportional to the prod- 
uct of the occupations in the two sublattices, namely 
4n d (l - n d ) ~ 1 - An 2 . 

It should be stressed that standard approaches that 
incorporate Coulomb interactions into the polaronic con- 
ductivity would predict an additive contribution VAn/2 
to the activation barrier at the CO transition^ which 
naturally arises from the extra energy cost that a po- 
laron in the majority sublattice has to overcome to hop 
to a neighboring unoccupied site. While such contribu- 
tion is indeed recovered in our treatment [see Eq. (|B13|) 
in Appendix B], we see here that the dominant conduc- 
tion mechanism at finite temperatures is related to the 
presence of charge defects. Being already thermally ex- 
cited onto the minority sublattice, such defects do not 
pay an extra cost in Coulomb energy when hopping from 
site to site, and their motion is therefore more advan- 
tageous than ordinary hopping across the gap (remind 
that at large A the number of defects n d does not depend 
explicitly on V, as was shown in Fig. [3]). 

Eq. ([M]) is compared with the full DMFT results for 
A = 1,V = 1 and A = 2,V = 1 in Fig. [5]b and Fig. 
IHc respectively. It should be stressed that the strong 
coupling analysis underlying Eqs. (|2"5j) and (j2"4"|) above is 
strictly valid only at A — > oo, as it neglects corrections 
proportional to D in the activation barrier. As a result, 
the above formulas widely underestimate the electrical 
conductivity^ and the agreement with the DMFT data 
at finite A and V — is not as good as the one obtained 
for cr(u>) in the preceding Section. This is particularly 
evident when comparing Eq. (|23[) with the A = 1 data 
of Fig. Ob, where large deviations arise especially at low 
temperatures T <C Ep. Increasing V shifts the ordering 
transition to higher temperatures, apparently improving 
the accuracy of the approximation. 

V. DISCUSSION AND CONCLUSIONS 

We have calculated the optical and electrical trans- 
port properties in a model describing charge ordering in 
systems with strong electron-lattice interactions. The 
DMFT treatment used here properly accounts for fluc- 
tuations of the lattice polarization. This allows to un- 
derstand the role played by charge defects in the electro- 
dynamic response of the polaronic ordered phase. The 



effects of a moderate intersite Coulomb repulsion term, 
which is often present in interacting narrow-band sys- 
tems, have also been incorporated. 

In the weak electron-phonon coupling regime, we have 
seen in Sec. II that lattice fluctuations play an impor- 
tant role as soon as A > 0, as they strongly reduce the 
critical temperature compared to the mean-field BCS-like 
expcctationj 23 i 30 Our results show that these fluctuations 
also manifest in the spectral properties, as they provide 
a disordered effective medium that scatters the electron 
motion, leading to an appreciable broadening of the sharp 
spectral features predicted by mean-field theory. 

The effects of lattice fluctuations are however much 
more dramatic in the strong coupling regime, where the 
charge ordering transition occurs in a polaronic phase 
with a preexisting charge localization. In that case, po- 
laron formation moves most of the spectral weight to a 
finite frequency absorption peak, whose position hardly 
changes at the ordering transition. Instead, the distinc- 
tive signature of charge ordering in the strong coupling 
limit is a marked enhancement of the polaronic absorp- 
tion peak, that is directly related to the suppression of 
charge defects, i.e. charges promoted to the "wrong" sub- 
lattice by the thermal fluctuations. In the polaronic or- 
dered phase, such defects also govern the transport prop- 
erties, as their motion dominates over the usual hopping 
mechanism across the disproportionation gap. The com- 
peting scenarios described above appear to coexist at in- 
termediate electron-lattice interactions (A > 1), i.e. pre- 
cisely in the region of the phase diagram where the crit- 
ical temperature is maximum. In this region, polaronic 
charge ordering reveals in the optical spectra through 
both a transfer of spectral weight to high frequencies and 
a global spectral weight enhancement as the temperature 
is lowered below T c . 

It should be noted that a model very similar to the one 
studied here was developed in a series of papers by Ihle 
and LorenZ ) 42 ' 44 i 45 and was applied to the long-standing 
problem of the Verwey transition in Fe304 — a proto- 
typical example of polaron ordering (see also the review 
paper Ref. |22j). In these works the model Eq. (1) was 
solved relying on the anti-adiabatic approximation for the 
lattice degrees of freedom, where the phonon dynamics 
is assumed to be faster than the electrons. Our adia- 
batic treatment can therefore be considered as comple- 
mentary to the one of Refs. I44ll45l . While the present 
approach clearly gives more insight into the physics of 
charge defects, our classical treatment of phonons has an 
important limitation that one should keep in mind when 
attempting a comparison with actual experiments. The 
results obtained here are expected to apply at tempera- 
tures and frequencies larger than the frequency scale of 
the phonons involved in the CO transition. Properly ac- 
counting for the quantum nature of the phonons would 
restore the possibility of a coherent transport regime at 
low temperatures^ 3 - and allow for the presence of charge 
defects even at zero temperature^ 
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APPENDIX A: MEAN-FIELD TREATMENT OF 
THE OPTICAL CONDUCTIVITY. 

We evaluate here the optical absorption in the CO state 
at mean-field level, which applies to the weak electron- 
lattice coupling limit. This treatment also describes ac- 
curately the numerical data at any finite V at nonzero 
temperature when A = 0, due to the intrinsic mean-field 
treatment of non-local interactions in DMFT. The mean- 
field solution for the electron propagators is obtained by 
inverting the matrix 



LO + A — £fe 

-e k w-A 



(Al) 



where A is the CO gap, which is determined self- 
consistently.— It is related to the charge disproportiona- 
tion by A = (XD + V/2)An. 

The corresponding spectral functions are delta func- 
tions, which allow for a direct evaluation of Eq. (|13|) for 
a generic DOS. The finite-frequency part, corresponding 
to interband transitions (across the gap) reads 



a(u) 



(2Af 



tanh ( — ) 
\4TJ 



2 v^ 2 ~ (2A) 2 

(A2) 

which has a square-root singularity at u> = 2A (the opti- 
cal gap) followed by a power-law decay at higher frequen- 
cies. At finite temperature, states are thermally excited 
across the gap, which enables the possibility of intraband 
absorption. In the absence of additional inelastic scatter- 
ing mechanisms, this gives rise to a zero-frequency peak: 



deN(e)(f>(e)- 



2c 2 



-/V 



RBZ 



-A 2 )] 
(A3) 



APPENDIX B: STRONG COUPLING 
APPROXIMATION TO o(uj) 

We derive here analytical approximations to Eqs. l[5]). 
([T4")) and (|2"Tj) that are valid in the strong coupling limit 
A> 1. 



1. Normal state 

Let us first note that the usual strong-coupling formula 
of Reikis for the optical confuctivity is recovered in the 



normal phase if we take the atomic limit (t = 0) in the 
propagators appearing in Eq. (|13[) . Indeed, by setting 



Gq~ = cj + i5 it can be shown that the lattice PDF 
is a sum of two gaussians centered at ±g/2k, each one 
carrying a weight ua — np — 1/3=2 (as was mentioned 
above, the Coulomb interaction term does not affect the 
properties of the normal phase at mean-field level). This 
can be cast in the form 



cosh(/?gA72) - 

r norm {JL)-g 



(gx) 



2 ^E 2 P 



AEpT 



(Bl) 



with the polaron energy defined as Ep = g 2 /2k = XD. 
The spectral function is obtained straightforwardly from 
Eq. dHU through the relation p no rm{u) = P n orm(u/g)/g 
[cf. Eq. ©]. 

Since in this limit the electron Green's function is site 
diagonal and momentum- independent, the e- integration 
can be factored out from Eq. (fT3")) . yielding a prefactor 
/ deN (e)(j)(e) = 1/4 on the Bethe lattice. The remaining 
frequency integral can be performed using the relation 



+ !/)_/(„)] 



sinh[^] 



2cosh[^±il)]cosh[^] 



(B2) 



which leads to the desired result 



a(w) 



<tq sinh(/3w/2) 



V / 2irE P T 16lu 



8E P T 



(B3) 



Focusing on the finite frequency part this can be 

further simplified to 



(Jo 



32oj v / 2ttE p T 



(uj - 2E P ) 2 
8E P T 



(B4) 



This formula represents a modified gaussian absorption 
band having its maximum at lu = 2Ep and a width 
\/2s = 2y/EpT, which reflects the thermal fluctuations 
of the phonon field that couples to the electron motion. 
Note that if the phonon quantum fluctuations are cor- 
rectly taken into account, the width of the absorption 
band does not shrink indefinitely at low temperatures, 
but rather saturates to a finite value s 2 = 2EpujqM: One 
can observe that the numerical factor 1/32 in Eq. (|B4|) 
implicitly includes the density dependent factor n(l — n) 
expected from theories of independent polarons. 



2. Charge ordered state 

In the charge ordered state (Art ^ 0), the two con- 
tributions to the lattice PDF are unbalanced by a factor 
which takes into account the different fillings tia ^ np on 
nonequivalent sublattices. As a result, Eq. (|B1|) becomes 

Pa AX) = P norm (X)[l ± (1 - 2n d ) tanh(/3 3 A/2)] (B5) 
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where the + and — signs correspond to the A and B sub- 
lattices respectively. The lattice PDF is again the sum 
of two gaussian peaks, centered at the minima ±X$ of 
the adiabatic potentials Va,Bi whose weights are now re- 
spectively n d and 1 — n d - At A = 2, this approximation 
compares very well with the full DMFT result illustrated 
in Fig. [21 That figure also shows that for moderate V < 1 
the shape of the gaussian peaks in the PDF is unaffected 
by the electrostatic repulsion, which only enters implic- 
itly through the self-consistent determination of n d . 

From Eq. (|B5j) we obtain the spectral function for the 
two sublattices: 

... / VAn\ ( VAn 
p A {u) = (1 - n d ) g+ Iuj-\ — I +n d 9- Ioj -\ — 

VAn\ ,„ , / VAn 



p B {u) = nd g+ w - 



+ (1 - rid) 9- ^ - 



where 



■. exp 



aEp \ 



2s 2 



(B6) 



and a = ±1. 

Taking the low temperature limit (T <C \/EpT, D) in 
Eq. (TT3]) and using the spectral functions above, it is 
possible to obtain an expression for the finite frequency 
conductivity (uj 3> T). At V — it takes the simple 
form: 



a(uj) = 2[1 - 2n d (l - n d )]a norm (uj) 



(B7) 



which directly relates the observed enhancement of the 
optical absorption in the CO phase, in agreement with 
the general arguments presented in Sec. IIII CI In the 
presence of intersite repulsion, the above formula gener- 
alizes to 

(uj - 2E P - VAn) 2 
-~ ^E~pT (\-n d ) 2 



a(uj) 



l6u}^2nE P T 



(uj - 2E P + VAn) 2 

8^ (n d f 



(B8) 



which clearly reduces to Eq. (|B7|) when V = 0. Both 
expressions compare well with the numerical data of Fig. 
[5] at A = 2. Moreover, the total optical spectral weight 
obtained by integrating Eq. (|B8|) down to a cutoff fre- 
quency uj = T correctly reproduces the non-monotonic 
temperature dependence of the DMFT results shown in 
Fig. [Sic. Note that at extremely large V, this expression 
predicts a double-peak structure, which is not observed 
in the numerical data at moderate V. 

The results presented in this Appendix, derived under 
the assumption of a vanishing bandwidth, are strictly 
valid as long as s ^ DM- From the asymptotic strong 
coupling expression T c = Z)/(16A) + V/4 we can estimate 



D 



(B9) 



We see that at V = the condition s > I) is never 
realized in the relevant region around T c , and some finite 
bandwidth corrections to the polaronic lineshapes can 
be expected^ Examination of Figs. [5] and [5] shows that 
the accuracy of the strong coupling formula Eq. (|B7[) 
improves when an explicit Coulomb term V is included. 



3. d.c. conductivity 

The d.c. conductivity is obtained by taking the limit 
uj — * in Eq. (I13p . Some care must be taken in inte- 
grating the gaussian spectral functions in the presence of 
the factor l/[4Tcosh(i//2T)] originating from the Fermi 
functions. For T ^ Ep one obtains 



<Jd.c.(T) = 



(TO 



WnEpT 



(VAnT 
SEpT x 



(BIO) 



2n d (\ - n d ) 



TlEp 

2T 



+ (l-n d ) 2 e 



2T 



n 2 ,e 



2T 



which reduces in the normal phase to the usual formula 
for the polaronic conductivity: 



a d.c. K 1 ) 



(TO 



-E P /2T 



32 v / 2lrEp~T 3 / 2 



(Bll) 



In the CO state at temperatures T < T c , the leading 
contribution to Eq. (|B10j) is 

<? d .c.{T) = in d (l - n d )e- W^ar™(T). (B12) 

which shows that the dominant conduction mechanism in 
the polaronic ordered phase involves thermally activated 
defects. The exponential term is close to 1 at moderate 
values of V such as the ones studied here, and can be 
dropped. This expression breaks down as T C T c , when 
the number of defects n d vanishes. In this case a more 
conventional result is recovered 



&d.c.{T) 



00 



I&kEpT 



(gj, + VAn/2)^ 
2E P T 



(B13) 



This formula predicts an additive contribution V An/2 to 
the activation barrier, corresponding to the extra energy 
cost that a polaron has to overcome when hopping to its 
neighboring sites 42 



APPENDIX C: F-SUM RULE IN THE 
BROKEN-SYMMETRY PHASE 

In this section we demonstrate the f-sum rule for the 
optical conductivity evaluated through the DMFT for- 
mula Eq. (|13|). generalizing the demonstration of Ref. 
|47| to the broken-symmetry phase. We first rewrite the 
integral of the optical conductivity in the form 



<j(uj) 



(TO 

2 



de2V(e)0(e)J(e) 



(CI) 



RBZ 



with 
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dv 



— oo J — oc 



dv ,ft")-m Tr 



a x ImG e {v')a x ImG e (v) 



(C2) 



This expression can be transformed to 



dvf{y)Tr 2a x ImG e (u)a x ReG e (i') 



Tr f°° 

— / dvf(v)Im (j x Gt(v)a x G t (y) 

7T ./_~, 



Tr f°° 

— / dvf(v)Im 

J-oo 



dG,{v) 
' de 



The latter equality is obtained using the fact that 
d(G7 1 G e )/(ie = d(l)/de — and observing that from 
Eq. (J9]) one has dG~ 1 /de — —a x . Going back to Eq. 
(|C1[) and integrating by parts yields 



<r{w) = <T / de^^- I dvj(y)2p* B {v). (C3) 

JRBZ ae 



With the present choice of the current vertex <fi(e) for the 
Bethe lattice one has ^j^- = —eN(e) so that the above 
expression reduces to 



K 



(C4) 



where K is the total kinetic energy defined as 



K= I deN(e)e / dvf(v)2pf B (v). (C5) 

IRBZ 



Introducing the following spectral density, p tk (y) — 
—Im[y/ZAZB — £fc] _1 /7r, with za and zb defined after 
Eq. @, and the corresponding momentum distribution 
function 



(C6) 



which is the appropriate generalization of the usual rife 
to the broken-symmetry phase, the total kinetic energy 
Eq. (|C5j) can be rewritten as: 



K = 2 I deN(e)e n e = 2 ^ 



(C7) 



keRBZ 



* On leave of absence at the ICMM-CSIC, Madrid, Spain. 
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